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Abstract. Within the framework of linearized Einstein field equations we compute the gravito- 
magnetic effects on a test particle orbiting a slowly rotating, spherical body with a rotating 
matter ring fixed to the equatorial plane. Our results show that the effect on the precession of 
particle orbits is increased by the presence of the ring. 
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, Analogy between classical electrodynamics and the linearized Einstein field equations 

rj • has been largely studied (Misner, Thorne & Wheeler 1973). This analogy was recently re- 

• , vised by Franklin & Baker (2007) in the frame of a test particle orbiting a slowly rotating 

' sphere. In this work we have introduced a massive rotating equatorial ring surrounding 

II the sphere in order to simulate the effect of an axisymmetric mass distribution. 

O The gravito-electric and gravito-magnetic potentials obey a Maxwell-like set of equa- 

^ , tions (Franklin & Baker (2007)) directly derived from linearized General Relativity (Soffel 

^ • 1989). For a sphere with Newtonian angular momentum is — ^MR'^ljs, the correspond- 

I— I, ing potentials are given by: 

^ : . . GM G£s sine , 

>. Vsir) = , As(r) = — (1) 

ly-^ ' r 2 r'^ 

. On the other hand, for the ring with Newtonian angular momentum £r = ma^ujR, we 

' have for the potentials: 

g: Mr) = -Gm|:^g;a^P,.(c„.9)-^ 

$-H ' Restricting the discussion to the mass-monopole and spin dipole, the gravito-electric 

. , and gravito-magnetic potentials for the composite system take the form: 



Vir) = -G^Ji±^ , 
r 

A(r) = -^(£5+^H)^^, (3) 
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Figure 1. Spinning sphere with radius R, angular frequency uis and equatorial mass current 

The Lagrangian for a test particle in the field of these potentials (neglecting terms 
from the Schwarzschild problem) reads: 



r 2c r 



(4) 

Therefore, the conserved quantities are: the z-component of the angular momentum 
and the energy of the particle. Both quantities constraint the initial conditions for the 
test particle. 
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From equation ([4]) we obtain the equations of motion for the test particle: 



f = r9'^ +r sin OijP' 
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In order to get the trajectory of the test particle we conducted numerical integrations 
of equations ([51). It is important to note that the term related with the matter current, 
originating in the gravito-magnetic field, is a perturbation of the gravito-electric field. In 
fact, the gravito-magnetic terms are of order {v jcf' , i.e., the effect is important for large 
timescales. Gravito-magnetic effects onto the system are "amplified" by using c — \. 
We constraint our simulations to planar and closed orbits setting 9 = 7r/2 and ^ = 0. 
3D simulations will be shown in a forthcoming paper (Castaneda et al. 2009). For the 
orbits in the plane, following the Bertrand's Theorem (Goldstein, PoUe & Safko 2000), 
wc impose an initial angular velocity for the test particle to obtain elliptic orbits. Initial 
conditions are therefore those corresponding to the perihelion, (rp), where = 0, r = 
and (p are obtained from the first integrals of motion given by equations ([5]). 

Numerical simulations were conducted using a Runge-Kutta method implemented in 
Matlab 7.0 and compared with a 3D code written in C. From Figure the ring's gravito- 
magnetic contribution to the orbital precession is obvious, which is the central issue of this 
paper. Our main goal has been to study a system of astrophysical interest; our results are 
a first step to understand the main features from gravito-magnetic contributions. Gravito- 
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Figure 2. Left panel: Orbits of tiie test particle in the x-y plane for the test particle with the 
ring at rest. Right panel: Sphere at rest. Parameters used in the simulation: c = 1, M = lO^^fcg, 
R = 2* lO'^m, m = M/4:, a = 1.5R, Vp = 1007?, ra = lOOOJ? y nip = 1kg. 

magnetic effects, such as The Lense-Thirring one, follow directly from our simulations, as 
will be shown in Castafieda et al. 2009. Finally, Figure ^ shows a family of trajectories 
for different ring velocities for a fixed sphere rotation velocity. The upper panel gives the 
effect of corotation between the sphere and the ring; the lower one shows the effect of 
counter-rotating bodies. The result is the precession of the orbit as was shown in a special 
case by Franklin & Baker (2007). However, we note that the effect of the equatorial mass 
current is to increase the precession velocity. In fact, as the velocity of the ring increases, 
the effect on the precession is more noticeable. These results motivate us to extend the 
present problem to further studies. 
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Figure 3. Orbits of the test particle in the x-y plane for o^s = 2 * 10 *s ^ and lor — kus with 

k = -100,-10,-1,1,10, 100. 
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